arXiv:1506.00228vl [math.PR] 31 May 2015 


1 


New Results On the Sum of Two Generalized 
Gaussian Random Variables 

Hamza Soury, Student Member, IEEE, and Mohamed-Slim Alouini, Eellow, IEEE 


Abstract —We propose in this paper a new method to compute 
the characteristic function (CF) of generalized Gaussian (GG) 
random variahle in terms of the Fox H function. The CF 
of the sum of two independent GG random variahles is then 
deduced. Based on this results, the probability density function 
(PDF) and the cumulative distribution function (CDF) of the 
sum distribution are obtained. These functions are expressed in 
terms of the bivariate Fox FI function. Next, the statistics of the 
distribution of the sum, such as the moments, the cumulant, and 
the kurtosis, are analyzed and computed. Due to the complexity 
of bivariate Fox H function, a solution to reduce such complexity 
is to approximate the sum of two Independent GG random 
variables by one GG random variable with suitable shape factor. 
The approximation method depends on the utility of the system 
so three methods of estimate the shape factor are studied and 
presented. 

Index Terms —Generalized Gaussian, sum of two random 
variables, characteristic function, kurtosis, moment, cumulant, 
PDF approximation. 

I. Introduction 

The generalized Gaussian (GG) signals catched the interest 
of researches in the recent years. The GG distribution (GGD) 
is widely used to model additive noise and interference, 
since several perturbations may indeed have a non-Gaussian 
behavior. It should be noted that the GGD family has a 
symmetric unimodal density characteristic with variable tail 
length that depends on the distribution parameter, the so called 
shape parameter a. 

In Gl, one of the non-Gaussian families used to model the 
noise is the GG noise (GGN). In fact, the GGN is deployed to 
study the detection for discrete time signals in non Gaussian 
noise. Moreover, in Cl-lll the noise is modeled by GGN 
and the performance analysis of the communication system 
perturbed by such noise has been studied and different metrics 
(such as probability of error, bit error rate, probability of false 
alarm etc.) have been evaluated. More recent works on the 
performance of communication systems perturbed by GGN 
can be found in na-Gi, where the authors derived a closed 
form of the symbol/bit error rates of several modulations 
perturbed by GGN over a generalized flat fading channel. 
Furthermore, in H the multi user interference in Ultra Wide 
Band was modeled as GGN. In addition, Fiorina introduced 
in ii a method to estimate the distribution parameter a using 
an estimation of the kurtosis of the GGD that depends only 
on the second and fourth moments. 

From a mathematical point of view, let A be a random vari¬ 
able following a GGD with parameter a > 0, mean E[A] = ji, 

H. Soury and M.-S. Alouini are with King Abdullah University of 
Science and Technology (KAUST), Thuwal, Makkah Province, Kingdom 
of Saudi Arabia. The authors are members of the KAUST Strategic Re¬ 
search Initiative (SRI) in Uncertainty Quantification. E-mail: {soury.hamza, 
slim.alouini} @kaust.edu.sa. 


and variance E[(A — p)^] = cr^, i.e. X ~ GGD{fi, a, a). The 
PDF of X, noted fx{x), is defined in ||5] Eq. 2] as 

fx{x) = exp (-A°|a; - /r|°) Vx € R, (1) 

2r(l/a) 

where A = ^ = ^ '^he normalizing coefficient 

with respect to the shaping parameter a. E denotes the 
expectation operator. 

Note that the GGD can model different type of distribution, 
for example the Gaussian distribution is obtained by setting a 
equal to 2, while for a = 1 one get the Laplacian distribution, 
and for a —oo we get the uniform distribution. 

The complementary CDF (CCDF) of standard GGD (zero 
mean and unit variance) was given in ||5l Eq. A.4] and defined 
as the generalized Q-function Qa(')- E is given by 

= —:2_r(l/a,A“x“) Vx > 0, (2) 

2r(l/a) 

where r(a,x) is the upper incomplete gamma function 
Eq. 6.5.3]. Note here that Qa(x) = 1 — Qa(—x) for negative 
argument. Once the CCDE of standard GGD defined, the CDE 
of X, noted as Fx (x) can be expressed as 

Ex(x) = 1 — Qa ^ 

Other statistics of the GGD are investigated in what follows. 

Erom the previous state of art, we notice that the majority of 
works focusing on the GGD. Actually, the works on the sum 
of two independent generalized Gaussian random variables 
(GGRV) are not as many as those related to a single GGRV. 
In fact, the most recent work on the properties of the sum of 
two independent GGRV is given in Eol, where Zhao et al. 
studied the properties of the PDE of the sum of two GGRV. 
They proved that such PDE has the same properties of the 
PDE of GGD (symmetry, convexity, monotonicity...) but they 
did not compute the PDE of the sum. On the other hand, 
they gave an approximation of the PDE of the sum of two 
identically independent distributed (i.i.d.) GGD. In this work, 
we are investigating the PDE of the sum of two independent 
not necessarily identically distributed GGRVs and focus on 
finding a closed form for such a PDE. 

A. Motivation and Applications 

As mentioned in m, we need sometimes to study the 
sum of independent GGRV (SGG), or a linear system of GG 
white noise such as AR(1) process driven by a GG process. 
Eurthermore, in seismic signal processing, the received seismic 
signal is modeled by a convolution of seismic signals, where 
seismic reflective coefficient can be modeled by GGD HI. 
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Thereby the distribution of the sum, that appears also as the 
convolution of the single distribution, is needed in this seismic 
model. Moreover, in communications, it is shown above that 
in some instances the noise and the multi user interference can 
be modeled as GG white noise ID, na. Therefore, the total 
perturbations at reception, defined as the sum of noise and 
interference, is modeled as the sum the GG signals. Many 
other applications can be found in the literature to motivate 
the study of sum distribution of GG signals. 

Per consequence, it is important to study the statistics and 
the density of the SGG distribution. The PDF was not derived 
before and so the CDF. Actually, the approach used in this 
work is based on the CF which was investigated in Ha for 
a > 1. However in our case we are studying the CF for any 
value of a using another approach of calculation based on 
the properties of the Fox H function (FHF) m. From the 
CF of one GGRV, we get the CF of the sum. Hence, the 
relation between the CF and the distribution densities leads us 
to investigate the PDF and the CDF of such distribution and 
its statistics (moments, cumulant, kurtosis...). 


Alternative expressions of cos(x) and e®”* in terms of the FHF 
are available in lfT4l Eq. (2.9.8) Sc Eq. (2.9.4)] 


cos(a;) = \/tt H|^’° 


0.2 




4 

-1 


= -H 


1,0 

0,1 




(7) 

( 8 ) 


Hence the integral identity defined in 114] Eq. (2.8.4)], solves 
the integral of the product of two EHE function over the 
positive real numbers. As a consequence the CE can be re¬ 
written as 


in (t) - 


trbo 

^0,2 


^ ..2 


I 

_ I 

(o,i),(ii)J 


H, 


1,0 


V'^ itiiTrl,! 

r(l/a) 


0,1 

2 


Ax 


cr"r(l/a) 


dx 


4r(3/a) 


■r 


\ Q ! ’ Q ! / 

(0,1), (il) 


(9) 


B. Organization of the Paper 

The remaining of this paper is organized as follows. In 
Section II, we investigate the statistics of the GGD, such as the 
CE, the moments, cumulant, and kurtosis. In Section III, the 
PDF and some statistics of the SGG distribution are presented. 
Next in Section IV, we analyze a method to approximate the 
PDF of the sum by a PDF of one GGD and the performance 
of such approximation is studied. Finally Section V concludes 
the paper. 


Note that the result proved in (|9]l is valid for all positive 
shape parameter a > 0. A previous demonstration was derived 
in 03 but only for a > 1, while the authors provided 
an expression of the CF of GGD in terms of the Fox- 
Wright generalized hypergeometric function which 

is a special case of the FHF |[T4l Eq. (2.9.29)]. The CE 
expression and derivation presented in this work present one 
of the contributions of the paper. 


H. Generalized Gaussian Distribution Statistics 


A. Characteristic Function 

Let a > 0, and X be random variable (RV) following a 
GGD{p, a, a) 

Theorem 1. The CF of X, E[e®‘^], is given by 


ipa{t) 


r(l/a) 4r(3/a) 


( 1 --,-) 

(0,1), (il) 


(4) 


where H!\[ ] is the Eox H function (EHE) |[T4l Eq. (1.1.1)], 

US). 

Proof: Let us start from the definition of the CE and the 
PDF of the GGD. In fact the CF can be expressed in integral 
form as 

= E[e"^^] = [ e''^^fxix)dx 

Jr 

= / e*‘“^exp(-A“|xr)dx. (5) 

2r(l/a) Jr 

Since \x\ is an even function, the integral in Q is the cosine 
transform of the exponential component 

A 

Pa{t) = / cos(te)exp(-(Aa;)“)dx. (6) 

r(i/a) Jo 


B. Moment Generating Function 

The moment generating function (MGF) can be directly 
concluded from the CF by the relation Ma{t) = (pai—it), 
so the MGF is obtained by 


M„(<) 


prl.l -q-^r(l/a) n 

r(l/a) 4r(3/a) 


( 1 --,-) 

V a ’ a 2 

(0,1), (ii) 


( 10 ) 

In some special cases, the MGF can be expressed in terms 
of elementary functions. For example for Gaussian case, i.e. 
a = 2, and using the special case of the FHF in m Eq. 
(2.9.4)], the MGF of Gaussian is 


M2{t) = 


2 


[ 2 

(0,1) J 




( 11 ) 


Another special case appears interesting is the Laplacian 
distribution (i.e. a = 1). Using the identity |[T4l Eq. (2.9.5)], 
the MGE of Laplacian distribution is given by 


Mi{t) = 

= e*'"H}4 

ptM 


( 0 , 2 ) 

(0,l),(i 1) 


- (7 t 

( 0 , 1 )' 

2 

( 0 , 1 ). 




( 12 ) 
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C. Moments and Ciimulant 


Without loss of generality, we are focusing our study to zero 
mean random variables (i.e. /r = 0). 

Due to the symmetry of the PDF of GGD, the odd moments 
vanish and the even moments obtained as follows 


m2n{X) 

'm2n+l{X) 


= E[X2"] 

= E[W2"+1] 


_^ 

r(3/a)" r(l/a) 


^ j2n r(l/a)" r(i 


= 0 


(13) 

Once the MGF and the moments are obtained, one can 
investigate the expression of the cumulant generating function 
(CGF) and the cumulant of the GGD. Actually the CGF, 
Kx{t) (or Ka{t)), is defined as 


Ka{t) = logMa(f) 


= log( 


r(l/a) 


) 


+ logH};2 


-cr^r(l/a) 2 

4r(3/a) 


(0,1), (ii) 


(14) 


By definition, the n-th cumulant, noted kn{X), is the n-th 
term in the Taylor series expansion of Ka{t) at 0. 


Theorem 2. The even cumulants of a zero mean GG random 
variable X can be expressed in terms of the even moments of 
X by 


k2n{X) = - ^ 

rni-\-2rn2-\- ■ ■ -TLirin^n 


(2n)!(mi H-+ m„ — 1)! 


mi!TO2!... rrinl 


n 

l<j<n 


r(3/a)^T(l/a)(2j)! J 

and the odd cumulants are equal to zero. 


(15) 


Proof: A cumulant kn{X) is the n-th derivative of the 
CGF evaluated at zero, ^ 

Since the CGF appears as the composite of two functions, we 
may use the Faa di Bruno’s formula ifT^ Eq. (2)] that computes 
the n-th derivative of composite functions 





n\{mi -I- • • • -f rrin — 1)! 
mi!m2!... 

m^{t) i 



(16) 


n 

the sum is over mi,m 2 ,...,m„ such that ^^jrrij = n. 

i=i 

Thereby evaluating ( [161 l at zero and replacing the moment by 
its expression, one get the final expression of the cumulant. ■ 
The cumulant of low order are easy to expressed 


k2{X) = a^ 
k4{X) = a^ 

keiX) = a^ 


/ r(l/a)r(5/a) \ 

V r(3/a)2 ; 

V r(3/a)3 r(3/a)2 J 

(17) 


Another statistics appear interesting to evaluate, in occurrence 
the kurtosis. The kurtosis, Kurt{X), is defined as the fourth 
cumulant divided by the square of the second cumulant of the 
distribution. In the GGD case the kurtosis is equal to 


Kurt{X) 


kyx) 

k2{xy 


r(l/a)r(5/a) 

r(3/a)2 


(18) 


One can easily check that (fTSl) confirms that the Gaussian 
kurtosis is equal to 0 and the Laplacian kurtosis is equal to 3. 

At this stage, the statistics of one GGRV are expressed 
in closed form. The next section considers the densities and 
statistics of the SGG distribution. 


III. Sum of Two Independent GG Random Variables 

As known the CF of the sum of two independent RV is the 
product of their CFs. Since the CF of the GGD is defined in 
the previous section, the CF of the sum can be easily obtained 
and so the densities by inverse Laplace transform of the CF. 
In fact, let X ~ GGD{p,i,ai,a) and Y ^ GGD{ii 2 ,(r 2 , P) 
two independent random variables following a GGD, and let 
Z = X + Y he their sum. It is clear that the first and second 
moment of Z are easy to find 

E[Z] = p, = pi + p2 

E[{Z-tiy]=a^=al+al (19) 


A. PDF and CDF of the Sum of Two GGRV 

The PDF of a random variable is known as the inverse 
Laplace transform of the CL. The CL of Z is given by 


Pz{t) = px{t)pY{t) 


r(i/a)r(i//3) 


Hi,2 


o'ir(i/'^) ,2 

4r(3/Q;) 


(0,1), (il) 


H 


1,2 


4r(3//3) 


( 1 -^ 4 ) 

(0,1), (ii) 


By applying the Laplace transform inverse to 
of Z is given by the following theorem. 


( 20 ) 

the PDL 


Theorem 3. The PDF of the sum of two independent GG 
random variable can be expressed in terms of the bivariate 
FHF iZZl/ in (EB- 

The LHL of two variables iflTll . also known as the Bivariate 
Lox H-function (BLHL) H’ T T ’[•,•] is a generalization of 
the LHL. Its MATLAB implementation is outlined in ifTSll . 

Proof: Let A = and B = The inverse 

Laplace transform of the CL ( l20l i 


4r(3//3) 

of Z gives the PDL 


fz{z) = 


27r 


2r(l/a)r(l//3) 


Hi,2 


A 


H 


1,2 


B P 


( 1 - 14 ) 

(0,1), (ii) 


( 1 --,-) 

(0,1), (il) 


(23) 
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olni/a) olTil/P) 

(i 1,1), (0,1,1) 

(1--,-) 

^ a ’ Q; 2 

1 

Tf3/a)[z-grn2,/(3){z-gf 


(0,1), (ii) 

-1 

1—1 


Fz{z) 


1 

2 


sign{z — g) 

2r(l/a)r(l//3) "" 


oini/a) oini/p) 

(^1,1),(1,1,1) 

(1-^.^) 

1 

ni/a){z-gmi/f3){z-gY 


(0,1), (ii) 

(0,1), (ii). 


( 22 ) 


The first two FHF are even functions, so the integral becomes 
a cosine transform of the product of these two FHF functions 


fz{z) 


1 

2 r(i/a)r(i/^) 



A f 


( 1 --,-) 

^ a' 

(0,1), (ii) 


H 


1,1 

1,2 


B 


(i-if) 

(0,1), (il) 


cos {t{g 


z)) dt 
(24) 


As seen before, the cosine has a representation in terms of the 
FHF (|2l). So we are facing an integral that involves the product 
of three FHFs over the positive real numbers. Such integral is 
solved in ifTTl Eq. (2.3)] and it is expressed in terms of the 
BFHF, which give us the final expression of the PDF of Z. ■ 

The CDF of Z is the primitive of fz that vanishes at (—oo). 
Back to (l24li . it appears that the CDF is expressed in term of an 
integral involving the product of two FHFs and sine function. 
The latter can be expressed in terms of the FHF for positive 
argument ITT] Eq. (2.9.7)]. Thereby, the CDE of the sum of 
two independent GGRV becomes the integral of the product 
of three EHEs which is evaluated in terms of the BEHE. 

Corollary S.l.The CDE, Fz{z), of the SGG distribution is 
given in (12^ . 

In (i22l i. sign{x) gives the sign of the real number x. The 
results in Theorem^ and Corollary 13.11 represent new results 
and they were not investigated before, which make them the 
essential contribution of this paper. 


B. Statistics of Z 

In the following analysis, the zero mean case is considered, 
while the non zero mean random variable can be obtained from 
the zero mean random variable by a simple shift Z = Z^ + g.. 

1) MGF: As mentioned before the MGE of Z can be 
obtained from the CE by the relation Mz{t) = (fz{—it) which 
gives the MGE of Z as 


Mz{t) 


r(i/a)r(i/^) 


-cr^r(l/a) 

4r(3/a) 


H 


1,1 

1,2 


4r(3//3) 


(0,1), (ii) 

(1-if) ■ 

(0,i),(ii). ■ 


(25) 


2) Moment: The moments of Z can be obtained from the 
binomial formula that describes the integer power of the sum 
of two numbers. Hence known the moment of X and Y O, it 
appears that the odd moments of Z, m 2 n+i{Z), vanish while 
the even moments are given by 


m2n{Z) 


arr(i)" ^/2„\ /a?r(i)r(|)y' 


2fc + 1 


2n — 2fc + 1 


(26) 


3) Cumulant and Kurtosis: Erom the MGE, it is easier to 
get the CGE by applying the logarithm to the MGE Mz(t). 
Thereby the CGE of Z is the sum of the CGE of X and the 
CGE of Y, Kz{t) = Kx{t)+KY{t). Moreover, the cumulant 
of Z is expressed also as the sum of the cumulant of X and 
the cumulant of Y. Note that the odd cumulant are equal to 
zero while the even ones are given by 


k2u{Z) = k2u{X) + k2n{Y) 


E 


(2n)!(mi 


- 1 )!, 


m\-\-2m2-\- • • -nmn—n 

r(|)^T(i)(2j)! j 


mi\m2 \ ■ ■ ■ TO„! 



n 

l<fc<n 


r(|)^r(i)( 2 fc)! 


(27) 


Once the cumulant expression is evaluated, the kurtosis of Z 
can be expressed, per definition, in terms of the fourth and 
second moments as 

ki{Z) ki{X) + ki{Y) 


Kurt{Z) = 


\2 ■ 


k2{ZY {k2{X) + k2{Y)Y 

Note that k 2 {X) + k 2 {Y) = af + cr^ = Thus a relation 
between the kurtosis of Z, X and Y appears as 


Kurt{Z) = ^Kurt(X) + 2XKurtiY). (28) 
The final expression of the kurtosis of Z is thus given by 


On 


Kurt[Z) = 


<^fr(4)r(|) a|r(i)r(f) 


t 2^2 


r(f) 


3',2 


r(f)^ 


- 3 (29) 
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IV. Approximation of the PDF of the sum of two 
GGRV 

The expression of the PDF of the sum of two independent 
GGRV (EB is quite high complex since it is expressed in 
terms of the BFHF Therefore, an approximation of the PDF 
is highly recommended to simplify the calculations and study, 
in simple way, the performance of systems in which the PDF 
of the sum is needed, like, for example, the evaluation of the 
symbol error rate (SER) of an M-phase shift keying (MPSK) 
over an GGN channel. Such analysis needs the PDF and the 
CDF of the SGG distribution. 

In this section we are investigating the approximation of the 
PDF of Z by the PDF of another GG random variable with 
shape factor 7 to be determined. In nni, it has been proved that 
the PDF of the sum cannot be a PDF of one GGRV. However 
the authors proved that both PDFs have the same properties 
(symmetric, convexity, monotonicity...). Furthermore, the PDF 
of the sum of two i.i.d. GGRV was approximated by the PDF 
of GGD. From that analysis, an approximation of the PDF of 
Z by the PDF of GGD is needed and worth pursuing. 

As shown in ([T]i, 3 parameters are needed to characterize a 
GGD, namely, the mean, the variance and the shape factor. The 
mean and the variance are given in (fT^ . Therefore, we need 
to find a method to get the shape factor 7 . In what follows, 
three approaches are presented. 


A. Kurtosis Approach 

The first method to estimate 7 is by using the kurtosis of the 
distributions. Since the kurtosis of the sum is already known, 
the shape factor can be obtained by equalizing both kurtosis. 
Thereby, we get the following equation to solve 


Kurt{Z^) = Kurt{Z) 


r(i)r(f) 

r(f)^ 


r(fr(f)2 + ’ 

(30) 


while Z^ ^ GGD{ii, a, 7 ) is the approximated RV of Z with 
parameter 7 . 



7 


Fig. 1: The curve of ( 1 ( 7 ) for positive values of 7 . 


Lets define the ratio between the variance of X and Y as 
(5 = -j, so the equation on 7 can be written in terms of 6 as 

(T9 


r(i)r(f) 

r(f )2 


.2r(i)r(|) , r(i)r(|) 




r(|)^ 


+ 


r(|)= 



By knowing a, jS, and S, (l30l) is written as h{j) = G, where 
h{'y) is a function on 7 , and G is a known positive constant. In 
Fig.m the function h is drawn versus 7 to analyze its behavior. 
Therefore, it appears that the function h{-) is a bijection. As 
such the equation h{'j) = G has only one solution in the 
positive real axis. Which mean that 7 exists and it is unique. 
The value of 7 is given as jKurt in Table |T] for some scenarios 
along with other values of 7 obtained from other approaches 
that we will discuss later on. 


B. Best Tail Approximation 

Another method to estimate 7 consist of taking the best 
choice of 7 that minimizes the square error of the tail. In 
other words 7 is chosen so the error between the exact PDF, 
fz{z) and the approximated PDF fj{z) at the tail is minimal. 
The tail is defined so z is above some level z > na 

pOO 

iTati = argmin / {f^{z) - fz{z)f dz, (32) 

T>0 Jna 

where n is chosen to define the desired region of the tail of 
the distribution. The minimization in (l32l i cannot be solved 
analytically by the available tools since it contains a shifted 
integral of a BFHF and FHF which is not known yet. A 
numerical evaluation of "frail is given in Table |T] for different 
values of 6. 

C. CDF Approximation 

This method is used to obtain the shape parameter that 
minimizes the error between the CDF of Z and the approx¬ 
imated CDF. Such approximation will give an asymptotic 
approximation of the complementary CDF (CCDF) which 
is needed in the computation of the probability of error. 
Mathematically, the shape parameter is given by 

pCO 

ICDF = argmin / {F^{z) - Fz{z)f dz. (33) 

t>o Jo 

In Table |T] some numerical values of jcdf are given and 
a comparison between three methods of shape parameter 
estimation is available too. 

An overview from Table U shows that the optimal value of 
"frail is near the value given by the kurtosis for any value 
of n. It is clear also that "frail approaches closely to "fKurt 
specially for n = 2 for all values of S. This analysis confirms 
the use of the kurtosis to approximate the PDF of the sum 
of two independent GGRV by another GGD to get a good 
tail approximation, this may also confirms that the kurtosis 
measure the heavy tail. Unlike this observation, the 7 obtained 
by minimizing the CDF error is a little bit far from outcomes 
of the kurtosis method. To conclude, these three methods can 
be used according to the situation we are facing. 
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TABLE I: Shape parameter for the approximated PDF using 
kurtosis, minimum CDF error, and minimum tail error for 
CTi = 1 and different values of (a, /3, S) 


(«./3,<5) 

'iKurt 

ICDF 

iTail 

n = 0 

1 

2 

3 

( 0 . 5 , 0 . 5 , 1 ) 

0.626 

0.467 

0.768 

0.673 

0.624 

0.642 

( 0 . 5 , 0 . 5 , 2 ) 

0.604 

0.492 

0.762 

0.656 

0.603 

0.584 

( 0 . 5 , 0 . 7 , 2 ) 

0.633 

0.501 

0.861 

0.741 

0.636 

0.834 

( 0 . 5 , 1 . 2 , 1 ) 

0.779 

0.602 

1.160 

1.053 

0.757 

1.165 

( 1 . 5 , 1 . 5 , 2 ) 

1.673 

1.373 

1.738 

1.702 

1.683 

1.664 

( 1 . 5 , 2 . 5 , 1 ) 

1.908 

1.391 

1.979 

1.959 

1.952 

1.887 

( 1 . 5 , 2 . 5 , 2 ) 

1.753 

1.443 

1.842 

1.799 

1.771 

1.741 

( 2 . 5 , 3 , 3 ) 

2.295 

1.941 

2.226 

2.261 

2.267 

2.335 



Fig. 2: Exact and approximated PDF of the sum of two GGRV, 
for /3 = 1.5, (5 = 2, cTi = 1, and two values of a. 


D. PDF and CDF Simulations 

The illustrations in this section are made for /3 = 1.5, ^ = 2, 
and cTi = 1. In Fig.|2l the PDF of the sum distribution is drawn 
for two values of a (0.5 and 2.5) and p takes two values to 
split the curves of both cases. The exact and simulated PDF of 
Z are drawn in the same figure among with the approximated 
PDF. The latter is computed using the kurtosis and the optimal 
CDF methods. It is clear that the exact PDF matchs perfectly 
the simulated PDF. Far from the mean, the approximated PDF 
appears close to the exact PDF and both methods have a good 
tail approximation. For a < 2, the kurtosis and the optimal 
CDF method are close to each other and match only the exact 
PDF at the tail with huge difference at the mean as mentioned 
in cni. However, for a > 2, the kurtosis method presents a 
good approximation of the PDF even around the mean, while 
the optimal CDF method represents a good approximation of 
the CDF as it will be seen later. We omit the optimal tail 
method here because it is close to the kurtosis method as 
shown in Table J] However one can draw it easily using the 
values available in Table |T] 

Our second illustration, highlighted in Fig. [3 consists of 
drawing the CDF of the sum for a = 2.5, /? = 1.5, 5 = 2, 
and (Ti = 1 using all three methods to approximate the PDF 
in linear scale. It is noticed that the results obtained from the 
optimal tail method (for n = 3) are very close to those issued 



Fig. 3: CDF of the sum of two GGRV using the Kurtosis, 
optimal tail, and optimal CDF approximations methods, for 
a = 2.5, P = 1.5, 5 = 2, and cti = 1 


from the kurtosis method. Another observation is that all the 
methods are close to each other and close to the CDF at the 
saturation region, i.e. Fz{z) ~ 1. This result is more detailed 
in the next figure which shows the complementary CDF in 
Log scale. 



Fig. 4: Complementary CDF of the sum of two GGRV for two 
values of a. 

In Fig. m the complementary CDF of the distribution of 
the sum is drawn for two values of a (2.5 and 0.5). For 
both cases, the approximated CCDF using the optimal CDF 
method matches the exact CCDF. However, as seen in Fig. |2l 
for a < 2, the CCDF obtained from the kurtosis and optimal 
tail methods is not too close to the exact CCDF. While, for 
a > 1, they are close to each other and asymptotically close 
to the exact CCDF. 

V. Conclusion 

The statistics of the distribution of the sum of two inde¬ 
pendent GGRV were derived in closed form in terms of the 
FHF and BFHF respectively, and the distribution of sum was 
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approximated by another GGD using three estimation methods 
to get the shape parameter, namely kurtosis, optimal tail and 
optimal CDF depends on the treated case. Finally, all the 
methods present a good estimation of the shape parameter 
especially at the tail region. 
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